Micellization in the presence of polyelectrolyte 
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We present a simple model to study micellization of amphiphiles condensed on a rodlike polyion. 
Q ' Although the mean field theory leads to a first order micellization transition for sufficiently strong 

C/3 ' hydrophobic interactions, the simulations show that no such thermodynamic phase transition exists. 

4-^ . Instead, the correlations between the condensed amphiphiles can result in a structure formation very 

^ ' similar to micelles. 

a 

o , 

I. INTRODUCTION 

Interaction between polyelectrolytes and ionic amphiphiles has attracted significant attention in both the molecular 
[ biology and the condensed matter physics communities p]-[r^. The driving motivation for this surge in interest is 
Q.^ ■ the possible application of polyelectrolyte-amphiphile complexation in gene therapy ||^. One of the major stumbling 
ly-j ' blocks in designing a successful gene therapy is the lack of a transfection mechanism by which a DNA strand can 
^f) be inserted into a cell |^ . The problem arises as a result of strong electrostatic repulsion between the DNA and the 
molecular membrane, both of which are negatively charged. This repulsion prevents a DNA segment from coming in 
! contact with a cell membrane, thus precluding any possibility of transfection. To overcome the electrostatic repulsion 
■ a number of protocols have been developed. The most explored ones rely on genetically modified viral vectors. A 
number of complications which can arise from the viral gene therapy have stimulated a development of non-viral 
methods One such method explores the association between a negatively charged DNA and cationic lipids or 
surfactants. A particular method which has attracted much attention relies on formation of lipoplexes |^-||. These 
. are complexes composed of a DNA strand and cationic lipid vesicles. Unfortunately the non- viral methods are also 
r-| prone to problems since, as is well known, the cationic surfactants and lipids are toxic to an organism An 
Q interesting question which then arises is what is the minimum concentration of cationic amphiphile needed to form a 
O lipoplex or a surfoplex? This question is well posed, since it has been known for some time that the polyelectrolyte- 



ionic amphiphile complexation occurs in a cooperative manner [|10| , [11 14-1^. It is, therefore, possible to identify 
the location of the cooperative condensation with the critical concentration needed to form complexes. What is the 
internal structure of such complexes? Are the condensed amphiphiles uniformly distributed along the DNA or do 
^ . they form micellar aggregates on the surface of a polyion? As a first attempt to study this difficult problem we shall 



appeal to a very simple model ml 18 
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II. THE MODEL 



We consider a rigid polyion modeled as a cylinder of length Zh and radius Bb inside a uniform medium of dielectric 
constant e. The total charge of a polyion, —Zq^ is uniformly distributed along the length of the cylinder so that each 
one of the Z monomers has charge —q. To further simplify the calculations, the longitudinal and angular degrees 
of freedom are discretized. The surface of the cylinder is subdivided into Z parallel rings of n sites each. The m 
condensed amphiphiles are restricted to move between the Zn ring sites, on the surface of the cylinder. 

Each amphiphile has a charged head group and a hydrocarbon tail. For generality we shall take the head group 
to have charge aq. The hardcore repulsion between the amphiphiles requires that each site is occupied by at most 
one amphiphile. For the specific case of DNA and dodecyltrimethyl amonium bromide (DoTAB) at the cooperative 
binding transition m « 0.8Z 

We define the occupation variables (Tij, with i = 1, . . . , Z and = 1, . . . , n, in such a way that cjij = 1 if a surfactant 
is attached at j'th ring position of the i'th monomer and <Ty = otherwise. Since the number of amphiphiles is fixed, 
the values of occupation variables obey the constraint 

Z n 

^ ^ (7y = m . 

i=l 3 = 1 

The Hamiltonian for this model has the following contributions: 
• Electrostatic surfactant-monomers interaction: 

z n 

= -a^yy ^ =, (1) 



where we have introduced a dimensionless charge density, the Manning parameter p9|-p]|. 
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FIG. 1. Interaction b etween a surf actant in the i'th ring with the j'th monomer. The distance between the surfactant and 
the monomer is r = by^ B? + {i ~ j)^. 

• Electrostatic surfactant-surfactant interaction: 



m2 



E E - 



(3) 
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FIG. 2. Distance between a surfactant located on the i'th ring in the fc'th position and the surfactant located on the j'th 
ring and the Tth position, r — by^ (i — j)^ + 4_R2 sin^ {n\k — l\/n). 



• Hydrophobic interactions between the hydrocarbon tails: 

Z n Z-1 n 



(4) 



1=1 3 = 1 



i=l j 



where ai^n+i — <^ii- The first term of Eq. (^) is due to the interactions between nearest neighbor amphiphiles on 
the same ring, while the second term is due the interactions between equivalent sites on consecutive rings. The 
hydrophobicity parameter x can be related to the size of the amphiphiles alkyl chain [l4| . 

The hydrophobic interactions between surfactant and water produce an effective attraction between the amphiphiles, 
forcing them to stick together. By doing so they expel the water molecules from their vicinity, lowering the overall free 
energy. The pairwise additive form of the hydrophobic interaction adopted in Eq. is clearly an over simplification. 
Nevertheless, we expect that this simple expression will help shed some light on the structure of the polyion-amphiphile 
complexes. 



III. MEAN FIELD THEORY 



To begin our study of the distribution of surfactants along the polyion we shall appeal to the mean field the- 
ory 1^,^. A note of caution, however, must be raised. While the mean field theory is expected to work very well for 
Coulombic long ranged interactions, it might not be so successful with the short ranged hydrophobic forces. This is in 
particular so since the problem is intrinsically one dimensional, and the fluctuations associated with the short ranged 
hydrophobic interactions are expected to be significant jl^ . With this note of caution in mind we shall proceed with 
the mean field study. 

The Gibbs-Bogoliubov inequality puts an upper bound on the total free energy, J- < J-q + (H — 7^0)0 = ^, where 
J^o is the free energy associated with the trial Hamiltonian Ho- To perform the calculations wc shall take Hq to be of 
a particularly simple one body form, 

Z n 

T-Co = qYY<l>ipcrip ■ (5) 

i p 

The partition function associated with Tio can be calculated straight forwardly 

Zo = ^ cxp I -f3q ^ 4'^p(Jip > = J]^ {1 + exp i-Mip)} , (6) 

{cr} [ i,p ) i,p 

and the free energy is 
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The average occupation of site p on the ring z, pip =< aip >o, is then 



^o = -^^ln{l + exp(-/3#,p)} . (7) 



or 

1 . ^l 



\ pip 

The free energy associated with Tio can be rewritten as 



^) • (9) 



After evaluating the average of {H — 7Yo}o with respect to TYq, the upper bound to the total free energy becomes 



to 



, + - 2 y^(j-j)2+4i?2sin2(^lfc-Z|) 

Z—ln Z n /l \ 

i j i—lj i j \ r J / 

To find the optimum upper bound, Eq. (|Tf|) must be minimized with respect to the average site occupation, leading 

P'l = YT^^iM^ ' (^^) 



where /i is a Lagrange multipHer introduced to enforce the constraint j Pij ~ ™ ^^'^ 

Pki{} - SikSji) 



Z ^ Z n 



p^,, = E , + E E ■ 



fc v/^' + - ; ^(i-fc)2+4i?2 sin2(f|j-/|) 

-x(Pjj-i + Pi,j+i) - X [(1 - + (1 - S,i)pi-i.j] . (13) 

Using the constraint above, fi can be evaluated leading to a self consistent equation for the average site occupation, 

TO 

^ m + j:kA^-pki)eP'^^'^-^'"^ ■ ' 

Eq. ( p^ ) can be solved numerically to find the equilibrium amphiphile distribution along the polyion. 

In Fig. ^ we show the average number amphiphiles per ring along the polyion for ^ = 4.17 and x = 1. We note that 
the amphiphiles are uniformly distributed along the polyion except at the ends of the macromolecule, where their 
density is strongly depleted. 



4 




50 100 150 200 

X 



FIG. 3. Comparison between the profiles predicted by the mean field theory and simulations results, for ^ — 4.17, X = 1) 
Z = 201, m = 160, and n = 5. 

For X ~ 1-35 a curious phenomenon occurs, as shown in Fig. ^. At this value of hydrophobicity the sites of 
the central ring become preferentially occupied by the amphiphiles. This corresponds to a micellization transition, 
in which the strong hydrophobic attraction between the surfactants overcomes entropy to produce a mesoscopic 
aggregate of amphiphiles. As x increases, a number of other peaks appear. Within the mean field theory we find that 
the micellization transition is of first order. The fact that a short ranged interaction produces a first order transition 
in a pseudo one dimensional system should leave us concerned. To check the existence of this transition we have 
carried out a set of Monte-Carlo simulations (MC). 



FIG. 4. The mean field theory prediction for x =1-3 and 1.4, ^ = 4.17, Z — 201, m = 160, n — 5. The peak is an artifact of 
the mean field approximation. 



IV. SIMULATIONS 



To simulate this model we use a standard Monte Carlo with particle-hole exchange, not restricted to nearest- 
neighbors pairs. The density profiles and the energy were measured after thermalization, results being both time and 
sample averaged. The simulations considered where done for DNA with ^ = 4.17, Z = 201, n — 160, to = 5 and 
various values of x- The averages were obtained using 100 samples. 

For X < 1-35 we find that the mean field theory is in excellent agreement with the MC. For x > 1.35, on the other 
hand, the simulations do not find any evidence of the micellization transition present in mean field (see Fig. H). The 
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energy is a smooth function of x, with no indication of the first order micellization transition, Fig. g. As expected, 
the short ranged hydrophobic interaction can not result in a phase transition in a one dimensional system. 
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FIG. 5. The Monte Carlo density profiles for x = and x = 5 for a system with ^ — 4.17, Z = 201, n = 160 and m = 5. For 
X = 5, the mean field theory predicts the existence of many peaks, that are absent in simulation. 

In spite of the absence of a true phase transition, the correlations between the condensed amphiphiles can lead to 
formation of structures along the polyion. To study these, we have constructed a histogram of amphiphile cluster 
sizes within the Monte Carlo simulation. Here the size of a cluster is defined by the number of amphiphiles per ring. 
Fig. 1^ shows that for amphiphiles with short alkyl tails (small hydrophobicity) the clusters are composed of only one 
amphiphile, with larger aggregates being highly improbable. With the increase in x "we find, however, that this is 
no longer the case and a significant fraction of amphiphiles belongs to the maximum sized cluster of n amphiphiles. 
Although this is not a thermodynamic transition, the change in behavior evident in Fig. |^ can be associated with the 
micellization. 




FIG. 6. The energy dependence on x calculated using Monte Carlo, ^ — 4.17, Z — 201, n = 160 and m = 5. 
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FIG. 7. Histogram of cluster sizes for Z=201 and n = 160, ^ = 4.17. 
V. CONCLUSIONS 

We have studied a simple model of micellization in the presence of polyoloctrolyte. It is found that the moan 
field theory predicts a first order micellization transition for the ionic amphiphiles condensed on a polyion. This 
thermodynamic transition is an artifact of the mean field approximation and is the result of neglect of fluctuations 
associated with the short ranged hydrophobic interactions. Monte Carlo simulations show that the moan field works 
very well for small hydrophobicities, but fails completely for strong short ranged interactions. Indeed the simulations 
do not find any evidence of a phase transition. Nevertheless, if the hydrophobic interaction are sufficiently strong 
they will lc;ad to significant correlations between the condensed amphiphiles, which can be interpreted as a micellar 
formation along the polyion chain. 
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